Complementary reaction analyses and the isospin mixing of the 

4 states in 16 0. 

K. Amos0 S. Karataglidis0 and Y. J. KirrQ 

School of Physics, The University of Melbourne, Victoria 3010, Australia 

(Dated: February 9, 2008) 

Abstract 

Data from the inelastic scattering of electrons, and of intermediate energy protons and pions 
leading to "stretched" configuration 4~ states near 19 MeV excitation in 16 O as well as from 
charge exchange (p,n) scattering to an isobaric analogue (4~) state in 16 F have been analyzed to 
ascertain the degree of isospin mixing contained within those states and of the amount of <i 5 / 2 — p^ 2 
particle- hole excitation strength they exhaust. The electron and proton scattering data have been 
analyzed using microscopic models of the structure and reactions, with details constrained by 
analyses of elastic scattering data. 

PACS numbers: 21.10.Hw,25.30.Dh,25.40.Ep,25.80.Ek 



*Electronic address: amos@physics.unimelb.edu.au 
^Electronic address: kara@physics.unimelb.edu.au 
■^Electronic address: yjkim@cheju.ac.kr 



1 



I. INTRODUCTION 



The nucleus is an unique environment in that the strong, electromagnetic, and weak 
interactions all contribute in ways that are manifest in its static and dynamic attributes. 
Furthermore, as nuclear systems display a rich array of properties, it is not surprising that 
diverse reaction studies, using a variety of probes and processes, are required to provide a 
range of complementary information before a detailed understanding of nuclear structure is 
possible. A requirement is that the reaction mechanisms of importance in each use (types 
and specifications) are well known. 

With the electromagnetic interaction, for both 7-decay probabilities and electron scatter- 
ing form factors that supposition is well founded at least for momentum transfers to ~ 
3 fm -1 . At higher momentum transfer values, meson exchange current (MEC) corrections 
may have noticeable effect particularly on transverse form factors. For electric transitions, 
other studies 0, |j| have shown how effective operators can be defined to account for MEC 
effects. With large basis shell model calculations of structure those corrections have led to 
good agreement with E2 transition data 0, 11 • Magnetic transverse form factors also are 
affected by MEC, but as their evaluations explicitly involve the nuclear currents, individual 
MEC diagrams must be evaluated. While the one pion exchange current (seagull and pion 
in flight) contributions have little effect at small linear momentum transfer values in scat- 
tering to states involving large angular momentum change, they have some effect at higher 
momentum transfer values and do so in the case of the M4 values in 16 p, though the 
effects are still minor and lie within vagaries of the choice of single nucleon bound state wave 
functions Q. 

Data from the inelastic pion scattering from nuclei are of particular interest, especially for 
pions with energies that ensure the A resonance dominance of the underlying pion-nucleon 
(ttN) interactions effecting the transitions. The interest stems largely from the distinctive 
relative transition strengths of 7r + and ir~ interactions with protons and neutrons. As a 
consequence, ir^' inelastic scattering cross sections should be sensitive to the isospin of 
(and so isospin mixing in) the nuclear states. But the specifics of the irN g-matrices in 
nuclear matter are not well known. 

In contrast, credible nucleon- nucleon (NN) ^-matrices in the nuclear medium have been 
formed for an energy of an incoming nucleon from 25 to 250 MeV. Using effective interac- 
tions that map to those ^-matrices, microscopic model (nonlocal) nucleon-nucleus optical 
potentials have been determined with which very good predictions of cross sections and spin 
observables from the elastic scattering of nucleons with energies in that range and from 
nuclei that span the mass table J5j . Such model evaluations of proton elastic scattering 
from 6 ' 8 He and from 208 Pb gave estimates for the neutron skin thickness of 

each nuclei. Specifically, calculations of nucleon elastic scattering were made using complex, 
spin-dependent, and nonlocal optical potentials formed by folding an effective two-nucleon 
(NN) interaction with nuclear state density matrices from a large basis space shell model 
calculation. That effective NN interaction accurately maps to the NN (^-matrices half off 
the energy shell, solutions of the Brueckner-Bethe-Goldstone (BBG) equations for nuclear 
matter found for diverse Fermi momenta to 1.6 fm -1 and built upon the BonnB NN po- 
tentials. The nonlocal optical potentials so formed were used also to specify the continuum 
waves in a Distorted Wave Approximation (DWA) study of inelastic scattering. The effec- 
tive NN interaction was also used as the transition operator in those DWA calculations as 
was a large basis space shell model spectroscopy for the nuclear states. Remarkably good 
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comparisons with data resulted, requiring no arbitrary renormalization or adjustments to 
bound state wave functions [5j|. The inelastic scattering results so found, cross correlated 
extremely well with the electron scattering form factor predictions of the selected model of 
nuclear structure. 

The study of transitions to the 4~ states in 16 O is particularly interesting. First, the 
dominant spectroscopy of the residual states should be particle-hole excitations (d 5 / 2 — P^} 2 ) 
upon the 16 ground state. If an essentially closed shell model for that ground state is 
considered, only two such states should exist; one each with isospin (T) of and 1. But 



there are four known 4~ states in the adopted spectrum in the vicinity of 19 MeV 
excitation; the three lowest having observable excitation strength from the ground with one 
or more of the scattering of electrons, pions and protons. Two of those are dominantly 
isoscalar. The third is classified as an isovector state whose analogue in 16 F has been 
seen in charge exchange (p, n) experiments 0, 0, Furthermore, electron scatterin 



(transverse) form factors for two of those states have been measured and analyzed 
7T^ inelastic scattering cross sections for energies in the A resonance region have been 
studied (ljj . and the cross section from inelastic scattering and cross section and analyzing 
power from charge exchange scattering of 135 MeV protons have been reported |l3l l30l|. 

In the next section, the spectroscopy of the 4~ states is discussed giving a particle-hole 
model for their excitation (from the ground) and also a prescription to study configuration 
and isospin mixing among them. Then in Sec. Illll a brief outline of details of the calculations 
of electron, pion, and proton inelastic scattering and of charge exchange scattering to the 
isobaric analogue state (IAS) state in 16 F is presented. The results of those calculations are 
then shown in Sec. |IV|and conclusions we draw follow in Sec. |V| 



II. THE SPECTROSCOPY OF THE 4 STATES 

Simple shell model calculations of 16 readily predict two 4r states in the spectrum. 
These are the isoscalar and isovector combinations of the "stretched" particle-hole excita- 
tions upon the ground state with a yy 2 proton (neutron) moving to the d 5 / 2 orbit. The best 
of such simple model calculations [16| predict excitation energies for those states of 19 MeV 
in accord with observation. The transition strength for excitation of such "stretched" states 
however, depends upon the fractional occupancies of the p 3 / 2 and d 5 / 2 orbits in the ground 
state. So one must have a good description for the ground state as well. The fractional 
occupancies can have a marked effect upon the normalization of cross-section peak values as 
was noted in a study of the excitation of the isoscalar /isovector pair of 6~ states in 28 Si [l7| 
wherein it was shown that using a large basis projected Hartree-Fock specification of the 
ground state could halve the strength of transition obtained from using a packed shell model 
description. 

Diminution of particle-hole excitation strength may be anticipated with any "stretched" 
states isospin pair. But in the case of 16 0, there are three strongly excited 4~ states in the 
region of 19 MeV excitation. Specifically they have excitation energies of 17.775, 18.977 



and 19.808 MeV and (dominant) isospin values of 0, 1, and respectively [111] . There is 
a fourth listed now |llj at an excitation of 20.5 MeV to which we do not have scattering 
data. Additional isoscalar negative parity states are not unexpected though as such are 
noted in the evaluated spectrum as well as being predicted by shell model calculations that 
allow 3p — 3h excitations |l6[. But various questions about the distribution of (d 5 / 2 — p 3 j 2 ) 
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strengths among these states ensue, as does the question of how much isospin mixing occurs. 

Particle transfer reactions, e.g. 17 0(d, t) and 17 0(d, 3 He) to these 4~ states give some 
useful indication of the (d 5 / 2 — P~>/ 2 ) strength distributions [18j. Assuming that the reactions 
occur by direct population of the lp — lh components of the residual states, the spectroscopic 
factors (C 2 S) in the pure isospin limit must satisfy sum rules, namely 

^(C 2 S) ( d)t ) = C 2 £( d) 3 He ) , (1) 

and 

(C 2 S)( dit )\(T=l) = -C^d.SHe) • (2) 

These relationships are satisfied approximately by experiment, and the extracted values 
for each transition suggest that the isovector (18.977) state contains 97% of the isovector 
particle- hole strength while the summed result of excitation of the (isoscalar) 17.775 and 
19.808 MeV states would account for 92% of that strength. But extracted spectroscopic 
factors are very model dependent. Such percentages should be considered indicative at best. 



A. A particle-hole model of the 4 excitations in 16 O 



Denoting the ground state of 16 by | J 77 ; T) = |0 + ; 0), a simple model for the 4 ;Tf states 



is 



|4-;7»= (p^d-H-jT, 



a \u) X fi (|.i) 



|0+;0) 



(3) 



(M 4 ,M T/=0 ) 

in which the particle-hole operator has been coupled in both angular momentum and isospin; 
the brackets surrounding the subscripts on the creation and annihilation operators indicate 
that full coupling. This operator is defined by 



(M 4 ,0) 



3 5 

- - m 3 - m 5 



4 - Mi 



1 1 

2 2 



a — a 



T f ) a 



A 



a, 



(4) 



The summation is taken over all component particle projection quantum numbers, and 
Mt = as we consider an iV = Z nucleus. Herein we do not consider any reverse amplitudes, 
i.e. [pz — dj 1 ). They will have very much smaller probability amplitudes, and none if the 

2 2 

d§/2 orbit in the ground state is vacant. The normalization of both of these states then 
is fixed by the fractional occupancies Oj of nucleons of either type in the orbits (j) in the 
ground state. The normalization, as developed in Appendix |XJ is 



Af 



1 — (75 > U 3 



(5) 



The same procedure allows a model of an IAS in 16 F though an additional isospin pro- 
jection changing operator need be included. 
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For inelastic scattering (or charge exchange leading to an IAS state in 16 F) involving 
scattering from a nucleon in orbit ji with isospin y and then leaving a nucleon in orbit j 2 
with isospin x, the one body density matrix elements (OBDME), 



4~; Tf 



a h* ® a ny 



+ :0 



(6) 



are required. Note that, in this case, only angular momentum is coupled in the specification 
of the operator. That is designated by the absence of brackets around the subscripts of 
the creation and annihilation operators and the use of <g> as coupling operator. Clearly the 
transfer quantum number is I = 4 and the combinations of proton and neutron spectroscopic 
amplitudes must effect an isospin change equal to Tf. 

For inelastic scattering to the 4~ states in 16 then, as developed in Appendix IB 11 
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while for the charge exchange (p, n) reaction to the IAS 4" 
must be a neutron to effect the charge transfer, y = | ■ 
Appendix IB 21 



state in 16 F, 
= —x. Then, 
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as the struck nucleon 
and as developed in 

(8) 



B. Configuration and isospin mixing 



We presume that the ground state of 16 has good isospin (T = 0), but configuration 
mixing will certainly reduce the p 3 / 2 occupancy from the closed shell value (of 4). A large 
basis shell model calculation (in a complete (0+2) huo space) has been made of the 16 O ground 
state and that gave values for (2j + l)aj respectively of 3.87 and 0.14 for the 0p 3 / 2 and 
0c?5/2 shells. To describe the 4~ states however, we consider a simple three component 
basis space. The three states of that basis the a pure particle-hole isovector states and two 
isoscalar states, 



|4";T = 1) : 

|4~;T = 0) : 
I (np — nh) 4~ 
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(9) 



With these it is possible to create states having both configuration and isospin mixing. In 
particular we suppose with this basis that there is an optimal prescription for the three 
states of interest in 16 O namely, 
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(10) 



We presume that the basis state (c| cannot be excited by a simple particle-hole operator 
acting upon the ground state, so that the scattering amplitudes for the excitation of each 
state, Ei, then are weighted sums 



M(Ei) = C (Ei)M {b \ + C^EJMw 



(11) 
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where the matrix elements A4/ x \ are transition amplitudes for the excitation of the basis 
states (a\ (T = 1) and (b\ (T = 0). However, it is also convenient to express the M.(Ei) in 
a form 

M(Ei) = Nx {cos(ei) M T =o + sin(ei) M T =i} 
M(E 2 ) = N 2 {cos(e 2 ) M T =i + sin(e 2 ) M T =o} 

M(E 3 ) = N 3 {cos(e 3 ) M T =o + sin(e 3 ) M T=l } , (12) 

where N = \/ Cq + Cf and cos(e) = Ct/N. The coupling angle for each state e;, and the 
scales Ni, are to be determined from data magnitudes of the measured cross sections. For 
evaluations however, individual proton and neutron amplitudes will be formed whence 

M{E X ) = Ni {[cos(ei) + sin(ei)] M n + [cos(ei) - sin(ei)] M u } 
M{E 2 ) = N 2 {[sin(e 2 ) - cos(e 2 )] M* + [sin(e 2 ) + cos(e 2 )] M u } 
M{E 3 ) = N 3 {[cos(e 3 ) + sin(e 3 )] M n + [cos(e 3 ) - sin(e 3 )] M v } . (13) 



III. DETAILS OF THE SCATTERING DATA ANALYSES 

A particle-hole structure model for three 4~ states in 16 has been used in the past @, 
0, H(| to analyze transverse magnetic M4 form factors extracted from inelastic electron 
scattering measurements 0, 0] . So also have been the differential cross sections from the 
inelastic scattering of 164 MeV pions [l^, leading to those same three states. But 
the inelastic scattering (of protons) to those states has been considered only with an old 
phenomenological approach (2(1 0] . We reconsider all of the data but now with the proton 
inelastic scattering studied using a ^-folding model of the scattering and with a large basis 
shell model wave function describing the ground state. Details of the methods used in our 
analysis are given next in brief. 



A. Pion scattering 



The cross sections for the inelastic scattering of pions 7r =b were evaluated using a con- 
ventional, phenomenological, distorted wave impulse approximation (DWIA). The distorted 
waves were obtained using the optical model potential 22l |. 



V^\r) = -Z^b 0p p p (r)-N^b 0n p n (r) + Zb lp \V / >Jr)Y\ ■■■ .Y/>„, <| V • /,„( r)V >> . (| |) 



where the density function form is 
( \ 2 1 



Z-2 



cxp 




(15) 



with parameter values listed in Table HJ The inelastic scattering amplitudes also require 
specification of the in-medium interaction of pions with bound nucleons as can be developed 
from a direct reaction scattering theory. But as the irN ^-matrices in nuclear matter are not 
well established and, as the pion energy should give A dominance, a simple contact form for 
the free 7rN interaction (the Kisslinger interaction) has been assumed as the nuclear state 
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TABLE I: The parameter values of the n optical potentials 





7T+ 


7T 


CLp — Cl n 


1.805 


1.805 


bo P 


-1.81 + 0.41* 


0.41 + 0.32 * 


bon 


0.41 + 0.32 * 


-1.81 + 0.40* 


bi p 


5.70 + 14.0 * 


1.53 + 4.69* 


bin 


1.53 + 4.69? 


5.70 + 14.0 * 



transition operator to effect the excitation of unnatural parity states. Specifically, we have 
used 

Kn = Gq[o-{kx k')] <J(F W - r n ) . (16) 

The strength was taken by a match to that for the 18.977 MeV state excitation assuming 
that it is a pure T = 1 state as the particle transfer data suggest is most probable. That such 
is adequate for the interaction will be confirmed by the results we find from simultaneous 
analysis of the scattering of ir + and of 7r~ leading to the 4~ state excitation. These details 
are gross simplifications of specifics of the scattering process but suffice for use in assaying 
the amount of isospin mixing there is in the 4~ states. In the DWIA calculations, the single 
particle states have been represented by harmonic oscillator wave functions. 



B. Electron scattering 

Electron scattering form factor calculations have been made using the usual point-particle, 
one-body current operators corrected for finite particle size and target recoil. Relativistic 
kinematics were used to specify the momentum transfer values. The single particle bound 
states used in our first analyses were harmonic oscillators for an oscillator length of 1.7 fm. 
Subsequently Woods-Saxon (WS) bound state wave functions were used, with which other 
studies |3| have found improved shapes of calculated form factors. Thus we have used WS 
functions for bound nucleon states with binding energies (of the dominant orbits and in 
MeV) being -38.0 (0s 1/2 ), -14.5 (pp 3 / 2 ), -10.5 (0pi/ 2 ), -4.28 (0d 5 / 2 ), -3.87 (ls 1/2 ), and - 
3.59 (0d 3 / 2 ). Higher shell (small occupancy) orbits were included and defined with binding 
energies between —1.5 and —1.0 MeV; their exact values are not very important for they 
have little impact on results. 

Analyses of longitudinal form factors from elastic scattering of electrons, as well as of 
cross sections and spin measurable from the elastic scattering of protons, when the nucleus 
is treated microscopically, require the OBDME for the target ground state. Frequently they 
are just the nucleon shell occupancies in the ground state. A complete (0 + 2)fko shell 
model calculation of 16 gave the the OBDME numbers (the dominant values) listed 
in Table |H] The first two columns are the shell occupancies while the cross-shell values 
are given in the last three columns. When weighted by the ground state OBDME listed in 
Table ITT| the WS functions gave the densities that are shown in Fig. Those densities have 
root mean square radii of 2.61, 2.65, and 2.63 fm for the proton, neutron, and total mass 
densities respectively; in good agreement with the value of 2.62 fm [23|, and like that, a bit 
smaller than the experimental value of 2.712 ± 0.015 fm. 
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TABLE II: The dominant shell occupancies and cross-shell terms for elastic scattering 



(nlj)i 



occupancy 



(nlj 



(nl 



j)2 



0*1/2 
0^3/2 
QPl/2 



Of/ 

Of/; 

Is 



5/2 
3/2 
1/2 



1.9993 
3.8706 
1.8942 
0.1413 
0.0678 
0.0107 



0,9 



1/2 



0^3/2 
0pi/2 



1*1/2 
1^3/2 
lPl/2 



0.0334 
0.2036 
0.0587 



0.2 




r(fm) 



FIG. 1: (Color online) Proton (large dash), neutron (small dash), and total (solid) mass densities 
for 16 as defined by the WS bound state functions. 

Electron scattering data (form factors) for scattering on 16 O have been measured; the 
elastic form factors by Sick and McCarthy [24| and the relevant M4 form factors by Hyde- 
Wright et al. p. While the elastic data has been measured to momentum transfers up 4 
fm _1 , the M4 data span an effective momentum transfer range from 0.8 to 2.6 fm^ 1 . Of the 
three states considered, form factors from excitation of the (dominantly isoscalar) 17.775 
MeV state and from the (dominantly isovector) 18.877 MeV state have been well resolved. 
Only statistical upper bounds are known for the form factor from excitation of the 19.808 
MeV state. 

Finally note that we have not included any MEC corrections in our calculations of the 
M4 form factors. They have been considered in the past [t| and found to be relatively minor 
in effect. Therefore MEC should not alter whatever general description of structure we can 
deduce from comparison of evaluated form factors with use of the same structures in the 
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analyses of other scattering data. 



C. Proton inelastic and charge exchange scattering 

Cross sections for inelastic proton scattering exciting the 4~ states in 16 0, and for the 
charge exchange reaction to the 4~ states in 16 F, have been evaluated using a fully micro- 
scopic DWA theory of the processes ||. The distorted waves are generated from optical 
potentials formed by folding an effective in-medium NN interaction set for each energy 
from the BBG g-matrices with the OBDME of the target states 0. In coordinate space 
that effective NN interaction is a mix of central, two-body spin-orbit and tensor forces all 
having form factors that are sums of Yukawa functions. Then with the Pauli principle taken 
into account, optical potentials from the folding are complex, nonlocal, and energy depen- 
dent. Such are formed and used in the DWBA98 program |25( to predict elastic scattering 
observables. That same program finds distorted wave functions from those potentials for 
use in DWA calculations of the inelastic scattering and charge exchange cross sections. The 
transition amplitudes for nucleon inelastic scattering from a nuclear target Jj — > Jf have 
the form [^] 



r = Tf f s j mu ' v {yt 

X { j\k o 0)\(* JfMf (l---A)\ A ge// (0,l)Ai{|xt +) (k,:0)> I^Wl- ■ -A))} ,(17) 

where Q sc is the scattering angle and Aqi is the antisymmetrization operator. Then a 
cofactor expansion of the target states, 

■■■A)) = -^=J2 I'PimW) I*Jm(1, ■■■A)) , (18) 



VC4 z 



allows expansion of the scattering amplitudes in the form of weighted two-nucleon elements 
since the terms dj m (l) \^jm(1, • ■ - A)) in Eq. UHlare independent of coordinate '1'. Thus 

T = ^(*w(lr--^)|a| ama (l)aii»ni(l)l*JiM 4 (l,---^)> 

x( x «(k o 0)|(^ 2m2 (l)| g e //(0,l) Ai{|xt +) (k,:0)> |^ imi (l))} 
£ (_)0l ~ mi) 7!TTT ^ IM ^ N \ J f M ^ Oi^m! -m 2 \I -N) S$£P 
x(x { j\k o 0) (^ 2m2 (i)| ge// (o,i) Ai {\xl + \m) 1^(1))} , (19) 



where reduction of the structure factor to OBDME for angular momentum transfer values 
/ follows that developed earlier for excitation from ground of the 4~ states. 

The effective interactions 5 , e //(0, 1) used in the folding to get the optical potentials have 
also been used as the transition operators effecting the excitations (of the 4~ states). As 
with the generation of the elastic scattering and so also the distorted wave functions for use 
in the DWA evaluations, antisymmetry of the projectile with the individual bound nucleons 
is treated exactly. The associated knock-out (exchange) amplitudes contribute importantly 
to the scattering cross section, both in magnitude and shape 0. 
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Thus only the structure details are left to be specified. They are the spectroscopic 
amplitudes as have been defined in Sec. |Hj and the single nucleon bound state wave functions. 
Those are exactly the same as we have used in analyzing electron scattering form factors 
and as specified in detail in the previous subsection. 

IV. RESULTS OF CALCULATIONS 

Analyses of the data from the scattering of pions leading to the set of 4~ states are 
considered first since they most strongly provide evidence of isospin mixing in those states. 
Then, in the second subsection, we present the results of elastic scattering of electrons and 
of protons from 16 O. Thereafter we give results of evaluations of the M4 form factors from 
electron scattering while in the last subsection we present results of excitations of the same 
4~ states by inelastic scattering and charge exchange scattering of protons. 

A. Pion scattering to the 4 states in ie O 

The isospin mixing in the three states of interest is indicated most clearly by considering 
the peak cross section values in pion scattering. Specifically we made calculations of 7r =1= 
peak height values and plotted ratios of actual data values over those theoretical estimates, 
irrespective of whether those peaks occur exactly at the same values of momentum transfer, 
and with whatever overall strengths Ni were required to match observed peak values. Our 
base scale choice for the pion-nucleon strength Q was that which best fit the actual data 
from the excitation of the 18.977 MeV state. The data from ir + scattering essentially has the 
same magnitude as that from n~ scattering llal and that near equality strongly suggested 
that the state was purely isovector in nature |lq . 

With that interaction strength we then evaluated the cross sections from excitation of 
the other (dominantly isoscalar) 4~ states, ascertaining the normalization scales required to 
fit the data and allowing variation in isospin mixing. The variation of the peak height ratios 
Cexp/cth so obtained are shown in Fig. |2]as functions of the sines of the coupling angles e». 
Such reflects the amount of isospin mixing one might expect for each of the three states. 
A variation of the peak cross section ratio found with the scattering of electrons exciting 
the 19.808 MeV state is also shown as it also emphasizes the isospin mixing expectation for 
that state. The variation with isospin mixing of the 7r + scattering cross section peak ratio 
is markedly different to that for it~ scattering in all three cases and those variations are 
quite sharp. Where the two projectile variations cross then indicates the amount of isospin 
mixing to be expected. With the 19.808 MeV excitation, the electron scattering data are 
only known as upper bounds and so using those (small) values gives a very sharp change 
when isospin mixing in that dominantly isoscalar state is used to calculate the relevant MA 
transverse electric form factor. That variation confirms well the suggested isospin admixing 
with sin(es) = —0.17. Thus we assumed, initially as least, that the isospin mixing in the 
three 4~ states of 16 to be sin(ei) = 0.2 for the 17.775 MeV state, sin(e 2 ) = for the 
18.977 MeV state, and sin(e 3 ) = -0.17 for the 19.808 MeV state. 

Holtkamp et al. quote a ratio of 7r ± summed yields for excitation of the 18.977 MeV 
state as 0.964 ± 0.08 so allowing some small amount of isospin mixing. A conservative 
estimate is shown by the width of the dotted box around the cross over point in the middle 
panel of Fig. 
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FIG. 2: (Color online) Peak height ratios of experiment to theory for ir^ scattering to the three 
4T states in 16 O as functions of isospin mixing. 



In Fig. 01 the differential cross sections measured ^ for the inelastic scattering of pions 
from 16 and leading to the 4~ states are compared with DWIA calculated results. The 
interaction operator strength was set by assuming the 18.977 MeV state to be described by 
the pure isovector basis state \a) of Eq. Q. That is consistent also with previous DWIA 
analyses |26[ of the same data as well as of cross sections from the excitation of "stretched" 
6~ states in 28 Si. The fits to the data from the other states required scales of 0.3 and 0.34, 
for the 17.775 and 19.808 MeV states respectively when the isospin admixing suggested by 
the peak height ratio study described above. That isospin mixing was particularly crucial 
since the relative peak height ratios of the 7r + and tc~ vary so markedly with the mixing. 
However, the shapes of the cross sections calculated for the dominantly isoscalar states do 
not match the data as well as they could. Indeed a previous analysis |26| used spin transition 
densities extracted from electron scattering which resulted in better cross sections shapes 
than we show. But such improvements are not central to this composite data study and, in 
any event, would require a better specification of the 7r ± iV (7-matrices including allowance 
of the T = I transition amplitudes; components found to be of some import, despite A 

dominance, in studies of ^ asymmetries in cross sections for excitation of the | + state in 
13 C [27j. Furthermore, any such refinement has not been made since there are also vagaries 
with the phenomenology involved in the calculations. 
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FIG. 3: (Color online) Differential cross sections for ir^ scattering to the three 4~ states in 16 
with fits requiring overall scale factors as given in the text. 



B. Electron and proton elastic scattering from O 



The charge density shown in Fig. ^ is tested by using it in evaluations of the longitudinal 
form factor for elastic electron scattering. The results of using harmonic oscillators (b = 1.7 
fm) and WS functions are compared with the data |24( in Fig. HJ The data (filled squares) 
show that the oscillator (dashed curve) is a good result to 2.5 fm -1 while the WS result (solid 
curve does as well to the second minimum in data at 3 fm -1 . We note that a much better 
fit to these data above 3 fm -1 has been found by Mihaila and Heisenberg ji^] by using 
an exp(S') coupled-cluster expansion approach. But to achieve those remarkable results 
necessitated use of an effective 50hw space, two-body currents, and Coulomb distortion; 
things not included in the simple approach we have taken since our purpose is to correlate 
major effects in scattering of different probes with analyses of data for which momentum 
transfer values are less than ~ 3 fm -1 . 

A further test of the propriety of the structure assumed for the ground state is to use 
that in analyses of the elastic scattering of protons from 16 O. Those we have made using 
the ^-folding approach that has proved quite successful in recent years. Details of the 
approach are given in the review as are details of the Melbourne effective NN interaction 
that is actually folded with the spectroscopy. 

We consider proton scattering at 200 MeV first since past studies of data from many 
nuclei gave great confidence that the effective NN interaction at that energy was good, 
and that the method of analysis was appropriate. In Fig. 200 MeV data H from 16 
(filled circles) are compared with the results (solid curves) found using the g-folding method 
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FIG. 4: (Color online) The longitudinal form factor from electron elastic scattering on 16 O. The 
data [2J] are compared with the results of calculations made using harmonic oscillator (b = 1.7 fm) 
wave functions (dashed curve) and with the set of WS functions (solid curve) defined in the text. 

with the appropriate (200 MeV) effective NN interaction, WS bound state functions, and 
(ground state) OBDME from the complete (0 + 2)hu shell model. These predictions are in 
very good agreement with that data and especially the good description of the analyzing 
power in particular suggests that the ground state spectroscopy we have chosen is realistic. 

Using that spectroscopy and simply changing the effective interaction to that appropriate 
for 135 MeV protons lead then to the results shown in Fig. El Therein, the elastic scattering 
cross sections of 135 MeV protons from 16 O as a ratio to Rutherford scattering are compared 
with data ji^. The WS folding result is displayed by the solid curve while the dashed curve 
portrays that found by using the oscillator functions instead. There are some differences 
between these results but the WS result is preferred since the cross section forward of 60° 
is the most significant as a test of the model. 

C. M4 form factors from inelastic electron scattering 

The electron scattering data for excitation of the 17.775 and 18.977 MeV states are well 
defined but only a weak upper limit is known for the M4 form factor from excitation of 
the 19.808 MeV state. As shown earlier, the latter quite strongly confirms the degree of 
isospin mixing in that state. First we consider the form factors and the separate proton 
and neutron components that arise when excitations of the pure isospin basis states are 
considered. The results displayed in Fig. [7| were generated using WS wave functions. The 
separate proton (tt) and neutron (u) results are shown by the long and short dashed lines 
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FIG. 5: (Color online) Differential Cross section (top) and analyzing power (bottom) for the elastic 
scattering of 200 MeV protons from 16 O. 




FIG. 6: (Color online) Ratio to Rutherford cross sections for the elastic scattering of 135 MeV 
protons from 16 O. The solid and dashed curves display the results found using WS and oscillator 
bound state wave functions respectively. 
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FIG. 7: (Color online) MA form factors from electron scattering to the pure isospin "stretched" 
states of Eq. ©. 

respectively. The dot-dashed line is the pure proton (and neutron) form factor one finds 
by adjusting the weight values to get identical proton and neutron amplitudes. Clearly the 
underlying amplitudes interfere constructively and destructively respectively to yield the 
isovector and isoscalar excitation results. By small weightings, magnitude equal proton and 
neutron amplitudes can be found and that defines the isospin mixing required to give nett 
vanishing form factors under conditions of destructive interference. That occurs essentially 
with the excitation of the 19.808 MeV state. Likewise finite results for the M4 form factors 
from excitation of the other two actual states can be found by such mixing. However, 
in those cases, there is also configuration mixing to be considered. By selecting the isospin 
admixtures as suggested from the pion scattering results, the M4 form factors from excitation 
of the 17.775 and 18.977 MeV states by electron scattering are those shown in Fig. |H| The 
dashed curves are those obtained using oscillator functions whilst those found on using WS 
functions are displayed by the solid curves. Both sets of calculations required the scalings 
(down) from the values found assuming that there was no spreading of particle-hole strength 
over other 4~ states in the spectrum. The results are in reasonable agreement with those 
made previously [f| in which oscillator wave functions with oscillator lengths of 1.58 fm 
were used to improve the comparison of calculated results with the data. As we have used a 
larger oscillator length, our form factors are more sharply peaked but, as evident, they are 
consistent with results found using the WS functions and with the elastic scattering data. 

Specifically isospin mixing has been taken into account with the analysis of the 
17.775 MeV excitation with the mixing angle, sin(ei) = 0.2, while the 18.977 MeV tran- 
sition was taken as a pure isovector excitation. But a small isoscalar admixture in this does 
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FIG. 8: (Color online) M4 form factors from electron scattering to the 18.977 MeV and to the 
17.775 MeV 4" states in 16 0. 

not vary the peak magnitude markedly. Thus the pure isovector characterization suggested 
by the pion scattering is neither abrogated nor totally supported by these results. However, 
they do identify an overall normalization for the state of N 2 = 0.7. In contrast, the calcu- 
lated form factor for the (dominantly isoscalar) 17.775 MeV transition varies noticeably in 
magnitude with isospin admixture and using the amount selected from the 7r =l= scattering 
data analyses requires a normalization N\ = 0.36 when harmonic oscillators are taken as 
the bound state functions. But the dominantly isoscalar transition result arises from an 
adjustment to the amount of destructive interference between the proton and neutron M4 
amplitudes. Such makes assessment of the configuration mixing effects difficult. We place 
more reliability in the proton inelastic scattering results that we discuss next, especially 
since the separate (bound) proton and neutron amplitudes constructively interfere to define 
the cross sections. 

D. Proton scattering leading to the 4 states in ie O 

In this subsection, we present the results of analyses of 135 MeV inelastic scattering 
data [30j, of 135 MeV (p,n) data (to the IAS 4~ state in 16 F) and of the analyzing power 
measured in that experiment ^3[. The latter is of particular interest in that spin observables 
usually are quite sensitive to details in the calculations. As we shall see, at 135 MeV, the 
charge exchange data track the inelastic proton scattering cross section quite well with but 
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FIG. 9: (Color online) Differential cross sections for the inelastic scattering of 135 MeV protons 
from 16 O exciting pure isospin "stretched" 4~ states. 

a scale factor of 2. That should be so if the state in 16 F is a true IAS of the 18.977 MeV 
isovector state in 16 0. 

All calculations have been made using a DWA approach with the distorted waves being 
those generated from evaluation of the elastic scattering of 135 MeV protons from 16 and 
the transition operator being the same effective interaction that we used in the ^-folding 
solution of that elastic scattering. 

First consider cross sections from the excitation of the pure isoscalar and isovector 
"stretched" 4~ states. The results of those DWA calculations are displayed in Fig El The 
isovector and isoscalar results are those labeled by T = 1 and T = in both panels of the 
figure. Harmonic oscillator bound state functions with b = 1.7 fm were used to get most of 
those results. Shown in the right hand panel are the pure proton (it) and pure neutron (y) 
excitation cross sections. They are noticeably different and the interferences between them 
to effect the pure isoscalar and pure isovector excitations result in cross sections that have 
similar magnitude but different peak positions. This is a quite different phenomenon than is 
evident with the electron scattering form factors for which the isoscalar and isovector form 
factors have so vastly different magnitudes. In the left panel, the harmonic oscillator results 
(solid curve for T = 1, dot-dashed curve for T = 0) are compared with those obtained using 
the WS bound state functions (long dash curve for T = 1, short dash curve for T = 0). 
Clearly the differences are very minor until larger momentum transfer values are involved. 

The proton scattering results are determined by the nature of the transition operator 
and the characteristics of it that most influence the results are revealed in Fig. The 
cross sections from excitation of the pure 4~ isoscalar (left) and isovector (right) states and 
the contributions from the central parts of the transition operator (dashed curves), from 
the two-nucleon spin-orbit (L ■ S) components (long dashed curves) and from the tensor 
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FIG. 10: (Color online) Differential cross sections for the inelastic scattering of 135 MeV pro- 
tons from 16 O exciting pure isospin "stretched" 4~ states and the contributions from different 
components of the transition operator. 



(S12) components (dot-dashed curves) are shown. Clearly the dominant component in these 
4~ excitations is that of the tensor force while the L • S terms are more important for the 
isoscalar than for the isovector transitions. The interference between the tensor and L • S 
contributions in the isoscalar transition is the cause of the shape difference of the total result 
from that for the isovector excitation. 

The complete results for the 135 MeV proton scattering cross sections to the actual 4~ 
states in 16 are shown in comparison with data in Fig. ^2 The isospin mixing as suggested 
by the analyses of the pion and electron scattering data have been used and the results scaled 
to fit; the scalings being identified as the configuration mixing weights. The solid lines are 
the results of DWA calculations made using the WS bound state functions while the dashed 
curve displays the result for the 18.977 MeV excitation found using harmonic oscillators. 
Clearly the shapes of the measured data are well reproduced with the distinguishing features 
of the dominantly isovector and isoscalar transitions being most evident. The scales, iVj, of 
Eq. (|T2*j) that were required to obtain these matches to the cross-section magnitudes were 
iVi = 0.42, N 2 = 0.62, and N 3 = 0.49. 

Finally we consider the analyzin g po wer results to compare with the data taken in a 135 
MeV charge exchange experiment k|- The results are presented in Fig. in the left 



panel the total result as found using the WS bound state functions is displayed by the solid 
curve while that found using the harmonic oscillators is depicted by the dashed curve. In 
the right panel, the results found using the separate central (dash), tensor (dot-dash), and 
two-nucleon spin-orbit (long dash) terms in the transition operator are presented. Note that 
these component results are normalized against the individual component cross sections they 
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FIG. 11: (Color online) Differential cross sections for the inelastic scattering of 135 MeV protons 
from 16 exciting the three 4~ states. The 17.775 and 19.808 MeV data and results have been 
shifted by 0.1 and 0.01 in scale for visualization. The charge exchange (IAS) (p, n) cross section 
data divided by 2 are shown by the open triangles. 

form (see the right panel of Fig. fTU|) and so their relative importance should be considered 
with comparison of the component against the total cross section. However, the tensor and 
spin-orbit amplitudes in particular are needed and their interference is crucial in finding the 
final total result that matches the data as well as it does. 



V. CONCLUSIONS 

The simultaneous analyses of electron, pion, and proton scattering data leading to the 
4r states in 16 ascertain that the states at 17.775 MeV and 19.808 MeV are dominantly 
isoscalar having isovector admixtures specified by 

iVisiri(ea) =0.42 sin(11.5°) , (20) 

and 

iV 3 sin(e 3 ) = 0.49 sin(-9.8°) , (21) 

while the 18.977 MeV state is virtually pure isovector in nature with N 2 = 0.62 being 
expected from the proton scattering data analyses. Thus we assess that the fractional 
exhaustion of the (0c?3 — 0p7 ',T) strengths exhausted in these transitions are 

2 2 

(Qd3~0pl 1 ;T = 0) = |0.42 cos(11.5°)| 2 + |0.49 cos(-9.8°)| 2 = 0.40 

2 2 

(Oda - Opl 1 ; T = 1) = |0.62| 2 + |0.42 sin(11.5°)| 2 + |0.49 sin(-9.8°)| 2 = 0.40 . (22) 

2 2 
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FIG. 12: (Color online) Analyzing power for the charge exchange reaction of 135 MeV protons 
from 16 exciting the 4" IAS state at 6.41 MeV. 



These values are less than one half of what has been assessed in the past from particle transfer 
data. But as the transfer data assessment is very much phenomenological model calculation 
dependent, particularly with what overall strength of interaction is chosen in zero-range 
DWBA calculations, and with a limited model spectroscopy which ascribed J 71 = 2~ to the 
state at 17.775 MeV excitation. As this is now known to be a 4~ state, and there is at least 
a fourth 4~ state at 20.5 MeV excitation, we believe our deduced numbers to be the more 
credible. 
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APPENDIX A: NORMALIZATION OF THE |4~;T/) STATES 
The normalization constant M of the 4~; Tf states follows from 
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The ground state expectation values then reduce as 
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where the fractional occupancies (of each orbit) are defined from 



0" 



a ]m'a' a jma 



o- 



In this we assume that every orbit j, m is equally occupied with 
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where rij X is the number of nucleons of type x that are contained in the shell j in the target. 
For a packed shell, as rij X = (2j + 1), then Uj X = 1. 
Thus 
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and then, with charge symmetry assumed, the sum over bound nucleon isospins a gives just 
a factor of 2 and, 



M = V2 
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(A6) 



APPENDIX B: SPECTROSCOPIC AMPLITUDES 



We consider the general matrix element, from which the spectroscopic amplitudes to be 
used in DWBA98 code evaluations can be derived, namely 



M {x , v) = U-\Tf 



a nx » a hy 



(i) 

(M 4 ) 



where now we have 



a' j2X ® a jiy 



(i) 

(N) 



E (-)*■ 



mi) 



Jij2 r m 1 -m 2 



+ ;0 , 



/ N ) Oj 2m2X djimiy 



while 



\A- ] T f )=M 



1 (4T/) 



a (5 1\ X a (3 1\ 
^ 2 ' 2 V 2 ' 2 J 



(M 4 ,0) 



|0 + ;0> , 



(Bl) 



(B2) 



(B3) 



21 



with normalization as given in Eq. (jA6|) . 

Notice that the transition operator is coupled only in angular momentum and not in 
isospin. That is so we find spectroscopic amplitudes for bound proton and bound neutron 
excitations separately as are required in evaluations of both electron form factors and cross 
sections initiated by nucleons. Those amplitudes, in fact, are the reduced matrix elements 
from Eq. (|B"T|) . 



M^, v) = -(0I0X 



4M, 



4-;T f 
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5l4 &NM f Sji%4; 



(B4) 



Thus, evaluation of the spectroscopic amplitudes for these transitions is made using the 



inverse of this, namely 



3<5/4<5jvM / -A / f( Xi2/ ). 



For inelastic scattering, the nucleon type does not change so the x = y in the above. 
However to consider charge exchange reactions, for which x = —y, we have to revise the 
specification of the final nuclear state. That transition is to the IAS which relates to the 
T = 1 state in 16 by the action of an extra isospin projection changing operator. 



1. Inelastic scattering 

With x = y, the general matrix element expands as 
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Then using the fractional occupancy representation of the expectation, Eq. (jA3|) . this reduces 
to 
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so that, for inelastic scattering, 
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2. Charge exchange to the IAS 4 in 16 F 



First we need to specify the form that the IAS takes. With T = y^, m a\ ,o im i being 



the operator that changes a neutron to a proton, the 4 state in 16 F that is an IAS to the 
particle-hole model isovector 4" state in 16 has the form 
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where P is the normalization. Then, the particle-hole description of the 4 states requires 
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the IAS state becomes 
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Normalization then defines V since, with Oj mx aj mx is the V nucleon number operator 



leading to a JX 
1 = (A-;IAS\A-;IAS) 



3 5 , 



m, — m c 



> (0 + ;0|[l-cr 



4-M^(-)(H^ m3 - m5 | 4 -M^ 



ot /iOb /_iat 103 1 

§"i 3 2 2 m 5 2 f»5-2 2 m3 2 



1 — (75 1 

2 2 



+ ;0 



= p 2 - 



1 — (75 1 

2 2 



1 — (75_1 

2 2 



(73 1 

2 2 



(B12) 



so that on inversion and with identical fractional occupancies for proton and neutron shells, 
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The spectroscopic amplitude follows then by using Eqs. (|B11|) and (jB13|) in forming 
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This is v^2 times the value found for inelastic scattering as given in Eq. ()B7)1 and so the 
charge exchange cross section should be twice that for inelastic scattering to the IAS. 
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